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Abstract 

It is shown that that the Volume free electron laser is a self-phase-locking system. 
Equations are derived that describe both the phase-locking process and the depen- 
dence of generation on the times of bunches' entering into the VFEL resonator in 
the superradiation regime and for long pulses. The derived equations are applicable 
to the description of the stated processes in gratings formed by relativistic BWOs. 



1 Introduction 



The research into pulse power amplification in microwave generators is gaining 
importance nowadays [HIS]- It has been shown that even for short microwave 
pulses, the discharge processes in the resonators of generators place limitations 
on the radiation power amplification potential. The possibilities of coherent 
combining of fields generated by several microwave generators are actively 
studied in this connection as the means of tackling these limitations. Particu- 
larly, phase- locking is a subject of vigorous study [2H5]- 

It has also been demonstrated that two separate super-radiant backward wave 
oscillators (BWO) connected to one and the same voltage (power) supply of 
subnanosecond rise time can form two coherent waves, which causes a fourfold 
increase in the radiation power |5]. 

Thus, these experiments confirm that using number of BWOs in this way, 
two (or more)-channel nanosecond relativistic microwave generators can be 
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developed, whose total power will be as high as W ^ N"^ {N is the number of 
BWOs [1]. 

Study of the application potential of volume free electron lasers (VFEL) for the 
development of relativistic microwave and optical generators with increased 
power is another promising line of investigation in this field [6]. 

Two(three)-dimensional distributed feedback (DFB) arising in VFEL resonators 
with two-three dimensional spatially periodic structures (now often called elec- 
tromagnetic or photonic crystals) enables producing coherent radiation from 
wide electron beams or several beams (see Figs. [1], [2]). A key feature of the 
VFEL is that as a result of diffraction, the signal is transferred from one point 
to another, thus linking the points of the beam and making them generate 
coherently (see Figs. [H [2]). 
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Figure 1. 
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Figure 2. 

This means that the VFEL is a self-phase-locking generator of radiation (phase- 
locking inside). When number of electron beams pass through electromag- 
netic crystals, the radiated power will increase in proportion to N"^ as a result 
of self-phase-locking of the beams due to two(three)-dimensional DFB formed 
in a VFEL resonator. When beams are generated, there is always a certain 
spread in the times when the produced electron bunches enter a VFEL res- 
onator, which results from the instability of the current generation in a diode 
gap. 

In the present paper it is shown that despite this spread in the entry times 
of the beams, a two(three)-dimensional feedback formed in a VFEL resonator 
gives rise to the self-phase-locking process . Equations are derived that de- 
scribe the process of radiation generation by several electron bunches and 
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enable studying the process of superradiation produced by the bunches, de- 
pending on the difference between the times of their entry into the resonator. 
It is also shown that similar processes of self-phase-locking, which lead to the 
power increase W ~ A^^ can also be observed when generation from electron 
bunches is excited in a system of several relativistic BWOs, whose corrugated 
waveguides form a single resonator by means of, e.g., slots (long bridge corru- 
gated waveguides or diffraction grating) made in a cylinder along the axis of 
the waveguides. 



2 Easing equations for the system with a photonic crystal (diffrac- 
tion grating) with changing parameters 



In the general case the equations, which describe lasing process, follow from 
the Maxwell equations: 



- IdD Att^ - IdH 

TotH = -— + —J, rotE = --, 

c at c c at 

divD = 47rp, ^ + divj=0, (1) 

— * — * — # 

here E and H are the electric and magnetic fields, j and p are the current 
and charge densities, the electromagnetic induction Di{f,t') = J eii{r,t — 
t')Ei{r,t')dt' and, therefore, Di{r,uj) = eii{r,u)Ei{r,ijj), the indices i,l = 
1, 2, 3 correspond to the axes x, y, z, respectively. 

The current and charge densities are respectively defined as: 

j{f,t) = eY,v4t)6{f-fo,{t)), p{r,t) = eY.6{f- fa{t)), (2) 



where e is the electron charge, Va is the velocity of the particle a {a numerates 
the beam particles), 

^ = ^ \E{Ut),t) + -lUt) X HiUt),t)] - ^iUt)Eirait),t))] , (3) 
dt mja I c \ 
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here 7q, = (1 — %)~2 is the Lorentz-factor, E{fa{t),t) {H{fa{t),t)) is the 
electric (magnetic) field at the point of location of particle a. 

It should be recalled that can also be written as 

dpa d'JaVa f ^ , 1 



m — - — = e 



dt dt [ c 



E{f^{t),t) + -[v^{t) X H{f^{t),t)] , (4) 



3 



where pa is the particle momentum. 

Let us recall here that the change in the particle energy through the its inter- 
action with electromagnetic fields is described by the equation 

mc —— = eVaE{ra{t),t). (5) 
dt 



Combining the equations in ([T]), we obtain: 

-A^ + V(VB) + -j^^— ^. (6) 



The dielectric permittivity tensor can be expressed as e{f) = 1 + x(r'), where 
x(r) is the dielectric susceptibility. When x ^ I, can be rewritten as: 

^^^'"''^ - ^S/ t')mt')dt' = 47r (^^^ + Vp(f,t)] .(7) 



When the grating is ideal x(r) = J2t Xr(r)^*'^'^' where r is the reciprocal lattice 
vector. 

Let the photonic crystal (diffraction grating) period be smoothly varied with 
distance, which is much greater then the diffraction grating (ptotonic crystal 
lattice) period. It is convenient in this case to present the susceptibility x{f) 
in the form, typical of the theory of X-ray diffraction in crystals with lattice 
distortion [6]: 

x(r) = Ee^*^^^^X.(r1, (8) 



where $r(0 = / T{f')dl', f{f') is the reciprocal lattice vector in the vicinity of 
the point r'. In contrast to the theory of X-rays diffraction, in the case under 
consideration Xr can also depend on r. Moreover, Xt depends on the volume of 
the lattice unit cell f2, which can be significantly varied for diffraction gratings 
(photonic crystals), as distinct from natural crystals. The volume of the unit 
cell Vt{f) depends on coordinate and, for example, for a cubic lattice it is 
determined as VL{r) = (f^d2{r)d3{r) ' where di are the lattice periods. If Xrir) 
does not depend on r, the expression ([8]) converts to that usually used for 
X-rays in crystals with lattice distortion. 

Recall here that for an ideal crystal without lattice distortions, the wave, which 
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propagates in the crystal can be presented as a superposition of plane waves: 



E{f,t) 



oo 

E 

T=0 



-rT — Ujt) 



(9) 



where k^- = k + t. 

Let us now use the fact that in the case under consideration the typical length 
for the change of the lattice parameters significantly exceeds the lattice pe- 
riod. Then the field inside the crystal with lattice distortion can be expressed 
similarly to (Q, but with Af depending on r* and t and changing noticeably 
at the distances much greater than the lattice period. 

Similarly, the wave vector should be considered as a slowly changing function 
of a coordinate. 

According to the above, let us find the solution of ([7]) in the form: 

E{f, t) = Re l^e^('^-(""^-^*)| , (10) 



where (pri^ = lo k{r)dl + $r(^^; where k{r) can be found as a solution of the 
dispersion equation in the vicinity of the point with the coordinate vector r, 
integration is made over the quasiclassical trajectory, which describes motion 
of the wavepacket in the crystal with lattice distortion. 

Now let us consider the case when all the waves participating in the diffraction 
process lie in a plane (coupled wave diffraction, multiple- wave diffraction), 
i.e., all the reciprocal lattice vectors r lie in one plane. Suppose the wave 
polarization vector is orthogonal to the plane of diffraction. 

Let us rewrite ffTOl) in the form 

E{f,t) = eE{f,t) = eRe S^Aie'^^'^'^-'''^ + A^e'^'^^^'^''''^ + ...} , (11) 



where 



01 (r) = / k^{r')dl (12) 







02 (r) 



A;i(f')rf/ + y T{f')dl. 





(13) 
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Then multiplying ([7]) by e, one can get: 

AE(f,t)-i^ J e(f,t-t0i?(r,O^t' = 47re[i^^ + Vp(f,t)] .(14) 



Applying the equality AE{f,t) = V(V-E) and using f|TT]) . we obtain 

/^(^^ei(<^i(0-<^t)) _ e^('^i("^-'^*)[2iV0iVAi + iVA?i(r)Ai - fci(r)Ai], (15) 

Therefore, substitution of the above expression into (fT^ gives the following 
system: 



_gHV'u';--; [2iA;i(f)VAi + iVki{r)Ai - kl{r}Ai 



lp«(</'l{^)-^i) 

1 9a;^£_T-(u;, r) 9^42 



^ Scj dt 



^ ^( 1 (9j(f,t) ^ A 

+ conjugated terms = 47re — — h \/p{r,t) , 

ot ) 

\^i{4>2{^-^t)^2ik2{T)\I A2 + zV4(f)A2 - A;2(F)A2 
+— £o(c^, r)A2 + %- — + —e^{uj, r)Ai 



" ' ' (9a; c?t c 

1 dbJ^eriyi, r) dAx 



duo dt 



1 d'jir.t) - 

+ conjugated terms = 47re — — h \/p{r,t) 

Xc^ at , 



(16) 



where vector /c2(^) = ki{f) + T, SQ^UJ^r) = l + xo(^), here the notation Xo(^) = 
Xr=o{'f') is used, er{co,r) = Xt(^)- Note here that for a numerical analysis of 
flTBl) . if Xo 0, it is convenient to take vector ki{r) in the form fci(r) = 



Let us multiply the first equation by e^^^'^^^'')^'^*) and the second by e~*'^''^2{i")-'^t)_ 
This procedure enables neglecting the conjugated terms, which appear fast 
oscillating (when averaging over the oscillation period they become zero). 

Considering the right-hand side of ( !T6|l . let us take into account that micro- 
scopic currents and densities are the sums of terms, containing delta-functions, 
therefore, the right-hand side can be rewritten as: 
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e 



c2 



Here is the time of entrance of particle a to the resonator, Tq, is the time of 
particle leaving the resonator, ^—functions in f|T71) indicate that for the time 
moments preceding ta and following Tq, the particle a does not contribute to 
the process. 

Upon averaging the system of equations over the oscillation period — , we can 
write: 



2iki{r)VAi + iVh{r)Ai - ki{r)Ai + — eo(w,^^i 



imue 



did dt 



t+- 



where the phase (pi[rait' ,ta,rao),t'] = (f>i{ra(t' ,ta,raQ) — ^^Jt' ■ 



2ik2{r)WA2 + iWk2{r)A2 - kl{r)A2 + ^eo(w, 

1 (9w%(a;,f) (9^2 . 1 (9wV(^^, r) Mi' 

= / e^7„(^05(r-f«(^^^,,r-;o)e-^^^[-(*''*-^^-)'*'l^^(^'-t„)^(t, 

where the phase '^2[ra{t', ta, rao), = 4>2{ra{t', ta, fao) - oof and Eq = I + Xo, 

When several {N number) electron beams move through a spatially periodic 
medium, the sum J2a oyer the particles can be represented as a sum of con- 
tributions coming from individual electron beams to the total current: 

E = E + E+-E, 



where N is the number of electron beams. 
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Using the definitions of 0^ (see (fT2|) . (fT3|) we can obtain tlie following rela- 
tionship for the phases ^ra- 



Equations ©-(jS]), describing particle motion in electromagnetic fields, and 
equations (fT8|l - (fT9|l for the fields are written using slowly changing amplitudes 
A and phases = 4>m — ^t- They give a closed, nonlinear set of equations 
that defines the amplitude Am and phases (as well as the change in the 
energy of particles interacting with the fields) and can be numerically analyzed 
using, say, the large-particle method. 

Because of random distribution of particles in the bunches incident on a res- 
onator (electromagnetic, photonic crystal), the times of particle entry into 
the resonator as well as the distribution of the entry point coordinates f\, of 
the bunch particles over the entire surface of the resonator are random. Each 
bunch also has a certain distribution of initial velocities This enables one 
to average (ITSl) and ( fT9i) over the distribution of the quantities to-r^, and 
Such averaging can be made by generalizing the averaging method developed 
for the case of one-dimensional generators like TWT, BWO, PEL to the case 
of a non-one- dimensional distributed feedback (DFB) The equations obtained 
as a result of such averaging in a stationary case when one beam moves in a 
VFEL resonator are given in 

We shall further consider self-phase-locking arising when photons are emitted 
by several electron beams in a spatially periodic VFEL resonator in the case 
of quasi- Cherenkov (diffraction) spontaneous radiation mechanism (recall here 
that this radiation mechanism underlies the operation of conventional one- 
dimensional TWTs and BWOs). 

For better understanding, let us suppose now that a strong magnetic field is 
applied for beam guiding through the generation area. Electron beams move 
along the direction of this field. Let us choose the direction of beam motion 
as the 2;-axis. We shall also consider the case when the period of the res- 
onator's diffraction grating changes along the direction of the ^;-axis. In this 
case, equations (flSll and (fT9!) can be presented in the form: 




(20) 



and 




(21) 



2iki{r)VAi + iVki{r)Ai - kl{r)Ai—eo{uj, r)Ai 
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+^ 



. 1 duj'^eo{uj,r) dAi to 



du! dt c 
2 gi{r±:Z,t), 



+ —e-r\^, r)A2 + 



dt 



(22) 



where 



t+- 



9i 



N i aN 



2ik2{r)VA2 + iVk2{r)A2 - kl{r)A2 + ^e^iuj, r)A2 

. 1 duj'^£Q{uj,r) dA2 , ^ ^ -1 du}'^er{oj,r) dAi 

du dt ^ ^ ^ du dt 

8Triue'd2 / \ 
= 2 92{r±,z,t), 



(23) 



where 



N f CIN 

Here 0i(^.(t', W, <)) = C^''''"''''''-"^^h^,{z')dz' and 02(^a(t', W, O) = 

Jq"^^ '^"^'^"'^^^ k2ziz')dz'. If the period of the diffraction grating is constant 
along the 2;-axis, then 0i = A;i2;Q,(t', t^^j^, w^^^)- The sign ((...)) means averaging 
over the distribution of beam particles over the transverse coordinate of the 
entry points, the entry time, and over the velocity distribution of the beams 
entering the resonator; 



1? 



U^ 



1 



, ki = kT=o, k2 = ki + T. 



Let p{f±) denote the particle density in the plane transverse relative to the 
particle velocity and h denote the number of particles traversing the inner 
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surface of the resonator per unit time. We also make use of the fact that 



5{z - z^,{t',t^,,um = S{t' - T{z,ta,,um/ 



where 

dz' 



dz if' f 



dt' 



'^{^T^aNT'^al,) ~ ^oiN '''' J 



Here r is the time required for the particle entering the interaction area z = 
at time ta at initial speed -u^^^ reach the point z and MAr(z, ta^, -u^^^) ^^e 
speed at point z for the particle whose initial speed at z = and time 

equaled m^^' '^hile ^"Af^^^iv) jg ^j^^ particle speed at time t if at 2; = and 

tapf its speed was Obviously, for the expression to be finite, this speed 
should not vanish; otherwise, such transformation of the 5-function becomes 
invalid (this is possible, for example, in the case when the beam's current 
exceeds the so-called limiting current and the virtual cathode is formed). 

We shall further assume that the particles are not retarded significantly during 
the interaction and write the quantity in the denominator of = '^al- 
a result, we have = Pz{ta), where Pz{ta) is the beam's time-dependent 
density distribution along the z-axis. 

Let us also suppose that the duration of the bunches injected into the resonator 
is larger than the period ^ of excited oscillations of the electromagnetic wave. 
Taking into account the relationship for r = t^^ + — , which follows in this 
case from the (5-function, the distribution density Pz(t^) in this case can be 
removed from the sign of integration over the entry time t^Ar point t' — 
and so we have 

Pz{ta) = Pz{t' )■ 



In a real situation, the distribution of the beam's longitudinal velocity is much 
less than the average longitudinal velocity of the bunch, so we can remove the 
density p^ from the sign of integration describing averaging over the velocity 
distribution. However, analyzing the phase dependence on the velocity dis- 
tribution of particles in a bunch, one should bear in mind that the velocity 
distribution can appreciably affect the phase. 

As a result, we can obtain the following expression for gn- 
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N 

t— ^ I '^'^ 

XUN J e-^[^"(^(-'*°''^'"^))---(*-'''°')l fiuP)duPdto. (24) 
t — ^ 

Here dto means integration over the initial times, /{u^n) is the velocity dis- 
tribution function in bunch N, and ujy is the average velocity of the A^-th 
bunch; the right-hand part of equation differs from zero at times t > from 
the initial moment defined as the instant of time when the first particle of the 
first bunch enters the resonator. 

The derived set of equations enables describing the process of radiation from 
several beams in a spatially periodic system (photonic crystal), including the 
case when the beams move opposite to one another. The geometry when the 
beams move opposite to one another can be used for beam diagnosing in 
the bunch-bunch collision region in coUiding-beam storage rings. For short 
bunches, this set of equations describes the phenomenon of super-radiance 
produced when several bunches of relativistic particles pass through a VFEL 
resonator. Particularly, it is possible to investigate radiation as a function of 
the difference between the times of electron bunches entry into the photonic 
crystal and as a function of the transverse distance between the bunches mov- 
ing in the electromagnetic (photonic) crystal. 

According to the equations derived here, the electromagnetic field induced in 
the crystal by different beams does not contain random phases r^j^ and any 
longer, and the total field, as a result, is a coherent sum of the induced fields, 
which means that the radiation power increases asW ^ E'^ = {J2n -^n)^ — N'^Ei 
It should be noted that when the parameters Xt grow, \xt\ ^ 1) the plane- 
wave expansion of the solution to Maxwell'equations, which is used in the 
dynamical diffraction of waves in crystals, requires that for accurate descrip- 
tion of the radiation generation process a larger number of waves should be 
considered. However, in this case one can expand the electromagnetic field 
into the analogue of the Wannier functions, which are used for describing the 
band structure of electrons in crystals in the case of tight binding. 

Let us consider the following example. Let a spatially periodic resonator be 
formed by axially corrugated cylindrical waveguides, we shall choose the di- 
rection of the waguides' axis as the z-a.xis. Coupled through the slots in their 
walls (long bridge corrugated waveguides or diffraction grating) , the waveguide 
form a single spatially periodic electrodynamical system. In the general case, 
we have a 2D periodic system in the (,t, y) plane, orthogonal to the 2;-axis. 
The beams move along the z-axis. Depending on the position of the cylinders 
in the transverse plane, square gratings or more complicated structures can 
be formed, e.g the cylinders can be arranged in a circle (as it occurs in the 
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magnetron) . 



To describe generation in this system, it is convenient start with the expansion 
Maxwell's equations in terms of the eigenfuctions Yn{r±) of this transverse 
grating (see a similar approach used for describing the motion of fast electrons 
in crystals in channeling regime (mode)) [7]. In this case the eigenfunction 
^n(ri) is a sum of the localized Wannier functions Wn 



where k is the reduced wave vector, n is the set of indices defining stationary 
wave functions , e.g., the wave function used for the formation of the structures 
periodic in the transverse plane, and Rm is the coordinate of the centre of the 
elementary cell m of the structure periodic in the transverse plane. 

As a result, for the analysis of generation of radiation we obtain one dimen- 
sional along the 2;-axis equations, where the excitation current is the total 
current I produced by the beams moving in the system. Let us average the 
current I over the electron entry times, the distribution of the initial velocities, 
and the distribution of electrons in the transverse plane (which in this case 
have the peaks in the regions where the electrons from each beam producing 
the current I move). 

As a result we obtain the equations similar in form to those used for the 
analysis of the generation process induced by one beam moving in a waveguide 
that is spatially periodic along the 2;-axis axis (formed by , e.g., a corrugated 
waveguide of a relativistic BWO.) Hence, we can conclude that the considered 
system, excited by several beams, generates common coherent radiation. Now, 
let us give a more detailed consideration of the case when the resonator is 
formed by the two elements of the grating. 

When the resonator period is formed by corrugation of the waveguide sur- 
face, we obtain a system consisting of two corrugated waveguides coupled 
through, say, a slot. For the BWO in the stationary case when two stationary 
electron beams move through circular waveguides, this system was analyzed 
neglecting the influence of the wave moving in the same direction [8]. Using 
numerical analysis, the authors of [8j showed that at certain parameters, a 
single-frequency oscillation mode is stet in the system, i.e., in fact, coherent 
summation of the amplitudes of the fields induced by two separate beams is 
possible. As follows from the above analysis, such coherent summation is also 
possible in the case non-stationary excitation of the system by two pulses of 
electron beams. 

Note here that the equations derived in this paper enable taking account of 
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the influence of a coherent wave on the generation process in a system of 
several BWOs. Moreover, according to [6], just in the range of parameters 
where the amphtudes of the incident and diffracted waves are comparable, in 
a two-three dimensional periodic system the increment of radiative instability 
increases sharply and the threshold for the generation start drops dramatically. 

When applied to this case, general equations for describing the excitation of 
two relativistic BWOs by two pulses of electron beams can be written in the 
form: 

Neglecting dispersion in considering the generation process in a system of two 
BWOs with a constant grating period (corrugation period), one can write 
these equations in the form: 



2iki, 



'■ dz 



+ 2i 



7? dt 



+ 



■ c- i-2 1,2 Aa 

-fcO fi-mn l^lz ^1 



iW^P Aa _|_ o^^l Ah _|_ uP_^2 Ab _ Y 



(25) 



'^^k2z"-§^+'^^r-^- 

+ c2 ^-r^2 + c2 Aafe^l + c2 X.ab^2 — '-mnjlr 



^'ft-lz 92 T^^*c2 9t ^[c2^0 ft'mn '^IzJ ^1 



(26) 



2zfc 



OA!; 



dz 

2 



+ 2i 



:2 at ^ 



^0 l^mn '^2z 



A\ 



I i£f /lb_|_i£lvl /l^-l-iid-v^ 4'' — V l' 
+ c2 ^T^2 + c2 Aab^l + c2 X.ab^2 ~ '-mnJ2Ti 



We have for the BWO mode 

ylf)(^ = 0) = 0, Af)(^ = L) = 0, 



where L is the resonator length. 

For the case when more than two waveguides are involved, say, number - 
N number of pairs of equation sets are required, instead of the two sets given 
above, and the terms describing the waves produced by other waveguides that 
are similar to Xafe^i xlb^l ( l25|) . ( |26|) should be added to each pair of 
equations. 
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It is worth noting that in a real case of arbitrary Xt, the coefficients appear- 
ing in these equations should be considered as phenomenological coefficients 
determined from the experiment on the passage of an electromagnetic wave 
through such structures. 

The derived system of equations enables describing the process of genera- 
tion excited by the combined pulses of electron bunches in a periodic system 
of coupled periodic waveguides (artificial electromagnetic crystal, VFEL res- 
onator), which forms self- phase-locking of coherent oscillations. The derived 
set of equations enabled studying the dependence of the radiation power on 
the difference between the times of the bunches'entry into the resonator of 
such a periodic system. This equation set is applicable to describing radia- 
tion produced by bunches with various duration, and consequently in the case 
of short bunches it allows one to describe the phenomenon of superradiation 
and phase-locking in the system of several relativistic BWOs coupled into the 
grating. 



3 Conclusion 

The above analysis shows that when a spatially periodic system of a VFEL res- 
onator is excited by several (N) pulsed electron beams that enter the resonator 
at different times, the two-(three)-dimensional distributed feedback, formed in 
the resonator, gives rise to self-phase-locking of the radiation process and co- 
herent collective oscillations, which result into a square in N increase of the 
radiation power W with growing number of beams: W ~ N^. Such self-phase 
locking occurring in a VFEL resonator makes it possible for us to consider 
the VFEL as the self-phase-locking system. Using this equation set, one can 
describe generation of generation by several beams in different modes: as su- 
perradiation from several electron beams and as radiation from long beams. 
These equations also make it possible to study the process of generation of 
superadiation as a function of the difference between the times of bunches' 
entry into the resonator. 



References 

[1] Benford J., Swegle J. A., Schamilogly E. High power microwaves, Taylor and 
Francis, 2007, 531 p. 

[2] Clayborne D. Taylor, D. V. Giri, High Power Microwave Systems and Effects, 
Taylor and Francis Group, 1994. 

[3] W. Woo et al, J.Appl. Phys. 65, 2 (1989) 0021-8979/89/020861-06. 



14 



[4] C.B. Wharton et a/,Proc. of the 8th Int. Conf. on High-Energy Power Particle 
Beams (BEAMS'90), Novosibirsk, USSR, 1990, pp. 1229-1240. 

[5] A. A. Elchaninov et al, Zh. Tekh. Fiziki 81, 1 (2011) pp. 125-130. 

[6] V.G. Baryshevsky, LANL e-print larX iv:1211.47 69[physics. optics]. 

[7] V.G. Baryshevsky High-Energy Nuclear Optics of Polarized Particles, World 
Scientific, Singapore, 2012. 

[8] V.A. Balakirev, A.O. Ostrovsky, Yu. V. Tkach, Pis'ma Zh. Tekh. Fiziki 16, 19 
(1990) pp. 8-12. 



15 



